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ABSTRACT 

We construct partially ionized hydrogen atmosphere models for magnetized neutron stars in radiative 
equilibrium with surface fields B = 10 12 — 5x 10 14 G and effective temperatures T e g ~ a fewx 10 5 — 10 6 K. 
These models are based on the latest equation of state and opacity results for magnetized, partially 
ionized hydrogen plasmas that take into account various magnetic and dense medium effects. The 
atmospheres directly determine the characteristics of thermal emission from isolated neutron stars. For 
the models with B = 10 12 — 10 13 G, the spectral features due to neutral atoms lie at extreme UV and 
very soft X-ray energy bands and therefore are difficult to observe. However, the continuum flux is 
also different from the fully ionized case, especially at lower energies. For the superstrong field models 
(B > 10 14 G), we show that the vacuum polarization effect not only suppresses the proton cyclotron 
line as shown previously, but also suppresses spectral features due to bound species; therefore spectral 
lines or features in thermal radiation are more difficult to observe when the neutron star magnetic field 
is > 10 14 G. 

Subject headings: magnetic fields - radiative transfer - stars: atmospheres - stars: magnetic fields - 
stars: neutron - X-rays: stars 



1. INTRODUCTION 

Thermal radiation from the surface of isolated neutron 
stars (NSs) can provide invaluable information on the 
physical properties and evolution of NSs. In the last few 
years, such radiation has been detected from radio pulsars 
(see Becker 2000; Pavlov, Zavlin, & Sanwal 2002), from 
old and young radio-quiet NSs (see Treves et al. 2000; 
Pavlov et al. 2002), and from soft gamma-ray repeaters 
(SGRs) and anomalous X-ray pulsars (AXPs) (see Hur- 
ley 2000; Mereghetti et al. 2002a), which form a poten- 
tially new class of NSs ( "magnetars" ) endowed with su- 
perstrong (B > 10 14 G) magnetic fields (see Thompson 
& Duncan 1996; Thompson 2001). The NS surface emis- 
sion is mediated by the thin atmospheric layer (with scale 
height ~ 0.1 — 10 cm and density ~ 0.1 — 10 3 g cm" 3 ) that 
covers the stellar surface. Therefore, to properly interpret 
the observations of NS surface emission and to provide ac- 
curate constraints on the physical properties of NSs, it is 
important to understand in detail the radiative properties 
of NS atmospheres in the presence of strong magnetic fields 
(see Pavlov et al. 1995; Ho & Lai 2001, 2003a; Zavlin & 
Pavlov 2002 for more detailed references on observations 
and on previous works of NS atmosphere modeling). 

This paper is the third in a series where we systemat- 
ically investigate the atmosphere and spectra of strongly 
magnetized NSs. In Ho & Lai (2001, hereafter HL01), we 



constructed self-consistent NS atmosphere models in ra- 
diative equilibrium with magnetic field B ~ 10 12 — 10 15 G 
and effective temperature T e g ~ 10 6 — 10 7 K and assum- 
ing the atmosphere is composed of fully ionized hydrogen 
or helium; we focused on the effect of the ion cyclotron 
resonance at E = E Bi = 0.63 Y o (S/10 14 G) keV, where 
Y e = Z/A is the electron fraction and Z and A are the 
atomic charge and atomic mass of the ion, respectively, 
and showed that the spectra can exhibit a broad feature 
due to the ion resonance (see also Zane et al. 2001). In 
Ho & Lai (2003a, hereafter HL03; see also Lai & Ho 
2002, 2003), we studied the effect of vacuum polariza- 
tion on the atmosphere structure and radiation spectra of 
strongly magnetized NSs. Polarization of the vacuum due 
to virtual e + e~ pairs becomes significant when B > Bq, 
where Bq = ra\c 3 jeh = 4.414 x 10 13 G is the criti- 
cal QED magnetic field strength. Vacuum polarization 
modifies the dielectric property of the medium and the 
polarization of photon modes (e.g., Adler 1971; Tsai & 
Erber 1975; Gnedin, Pavlov, & Shibanov 1978; Heyl k 
Hernquist 1997), thereby altering the radiative scattering 
and absorption opacities (e.g., Meszaros & Ventura 1979; 
Pavlov & Gnedin 1984; see Meszaros 1992 for review). 
We showed that vacuum polarization leads to a broad de- 
pression in the high-energy {E > a few keV) radiation 
flux from the atmospheres and hence softer high energy 
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tails, as compared to models without vacuum polarization. 
The depression in the spectrum is broad because of the 
large density gradient in the atmosphere. This depression 
of continuum flux also strongly suppresses the equivalent 
width of the ion cyclotron line and makes the line more 
difficult to observe (for B > Bq). HL03 suggests that the 
absence of lines in the observed thermal spectra of several 
AXPs may be an indication of the vacuum polarization 
effect at work in these systems. 

In our previous studies, as well as in other previous stud- 
ies of magnetized hydrogen and helium NS atmospheres 
(see Pavlov et al. 1995; Zavlin & Pavlov 2002), the atmo- 
spheres are assumed to be completely ionized. Because 
the strong magnetic field significantly increases the bind- 
ing energies of atoms, molecules, and other bound states 
(see Lai 2001 for a review), these bound states may have 
abundances appreciable enough to contribute to the opac- 
ity in the atmosphere (see Lai & Salpeter 1997; Potekhin, 
Chabrier, & Shibanov 1999). Although the energy levels 
and radiative transitions of a stationary H atom in strong 
magnetic fields have long been understood (e.g., Ruder et 
al. 1994), complications arise from the non-trivial effect 
of the center-of-mass motion of the atom on its internal 
structure (e.g., Potekhin 1998, and references therein). 
Also, because of the relatively high atmosphere density 
(p > 1 g cm -3 ), a proper treatment of pressure ioniza- 
tion is important (e.g., Lai & Salpeter 1995; Potekhin et 
al. 1999). Recently, thcrmodynamically consistent equa- 
tion of state (EOS) and opacities for a magnetized, par- 
tially ionized H plasma have been obtained by Potekhin 
& Chabrier (2003, hereafter PC). In this paper, we study 
NS atmosphere models based on the EOS and opacities of 
PC (see Ho et al. 2003 for an initial report of the results 
shown here). 

Section 2 briefly summarizes the energy structure of hy- 
drogen atoms in strong magnetic fields. Some details of the 
atmosphere models are discussed in Section 3. We present 
atmosphere models and spectra for different magnetic field 
strengths and temperatures in Section 4. In Section 5, we 
examine a dense plasma effect not accounted for in pre- 
vious modeling of NS atmospheres. Section 6 summarizes 
and discusses the implications of our results. 



2. HYDROGEN ATOM IN STRONG MAGNETIC FIELDS 

In a strong magnetic field with b = B / Bq 1, where 
B = m 2 c e 3 c/h 3 = 2.3505 x 10 9 G, the magnetic field con- 
fines the electron gyromotion in the transverse direction to 
a size ~ do / In b, where ao is the Bohr radius, and the elec- 
tron is restricted to the ground Landau level. The hydro- 
gen atom energy spectrum is then specified by two quan- 
tum numbers (s,v), where s measures the mean trans- 
verse separation between the electron and proton (or the 
guiding center distance of the electron gyromotion) and 
v gives the number of nodes in the z wave function. For 
ln[b/(2s + 1)] ^> 1, the energy of the states with v = is 
Eg = — yt s (ln6) 2 , where the superscript "0" indicates the 
energy of an atom with a non-moving (infinitely massive) 
nucleus and A. s depends weakly on b and s and formally 
tends to 1 Ry (= 13.6 eV) at b — > oo. For example, for 
s = and 10 11 G < B < 10 16 G, A. = (5.0 ± 0.7) eV. 
A convenient fitting formula for arbitrary B (accurate to 



within a few parts in 10 3 for B < 10 17 G and s < 3) is 

E° s = -A s x 2 , x = ln(l + cob), (1) 
A a _ 1 + c 2 x- 2 + c 3 x~ 3 + c 5 x~ 5 + c 6 (l + s)~ 2 x~ 6 
Ry ~ 1 + ax' 1 + c 4 x~ 2 + c e x- 4 ' ^ ' 

Here, c 5 = c 6 (l + s)/co ensures the correct behavior of the 
fit at b <C 1, and the parameters a (i = 0, 1, 2, 3, 4, 6) are 
given in Table 1 for s = 0, 1, 2, 3. 

For the states with v > 0, the energies at fc > 1 are 
Eg U ~ n~ 2 Ry, where n is the integer part of (y + l)/2 
(see Lai 2001, and references therein; see also Potekhin 
1998 for a more comprehensive set of fitting formulae). 

This simple picture of the H energy levels is modified 
when a finite proton mass is taken into account. Even for 
a "stationary" H atom, the energy E s becomes E®+sHu)b p , 
where tv-OBp — 6.3 B\2 eV is the proton cyclotron energy 
and B12 = B / (10 12 G); the extra "proton recoil" en- 
ergy sfuvBp becomes increasingly important with increas- 
ing B. Moreover, the effect of center-of-mass motion is 
non-trivial: When the atom moves perpendicular to the 
magnetic field, a strong electric field is induced in its rest 
frame and can significantly change the atomic structure. 
Indeed, in the limit of large transverse pseudomomentum 
p±, the atom assumes a decentered configuration, and its 
energy depends on p± as — pj 1 rather than the usual p\ 
dependence (see Potekhin 1998 for fitting formulae and 
references). 

The most important transitions for the radiative opac- 
ity are the bound-free transition from the (s,u) = (0,0) 
state and the bound-bound transition (0,0) — > (1,0). The 
former produces an absorption edge at the energy E$, 
which is smoothed out due to the dependence of Eq on p± 
from magnetic broadening (see Potekhin & Pavlov 1997) 
and due to the merging with transitions to v > states 
at slightly lower energies (PC). The (0,0) -> (1,0) tran- 
sition produces the maximum absorption at the energy 
E (p± - 0) - £i( P _l = 0) = Eg - E° + hjj Bp - with in- 
creasing Eo{p±) — E\{jp±) first decreases to a mini- 
mum and then increases to fiws p . As a result, this bound- 
bound transition produces a magnetically broadened fea- 
ture that resembles an inverted absorption edge, rather 
than a line, plus peaks at min(.Eo — E\) and huJBp (Pavlov 
& Potekhin 1995). The transitions from s = to s > 1 
are forbidden for stationary atoms but become allowed for 
moving atoms and thus contribute to the opacity. See PC 
for more detailed discussion of opacity features. 

3. MODEL METHOD 

All models considered in this paper have the magnetic 
field aligned perpendicular to the stellar surface, and thus 
the spectra presented are emission from a local patch of 
the NS surface. We solve the full, angle-dependent radia- 
tive transfer equations for the two coupled photon modes 
in order to construct self-consistent NS atmosphere mod- 
els. The two photon polarization modes in the magne- 
tized plasma that characterizes NS atmospheres are the 
extraordinary mode (X-mode), which is mostly polarized 
perpendicular to the k — B plane, and the ordinary mode 
(O-mode), which is mostly polarized parallel to the k — B 
plane, where k specifies the direction of photon propa- 
gation and B is the direction of the external magnetic 
field. In our models, the temperature corrections AT(r) 
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Table 1 

Parameters for energy levels of hydrogen atom [see eqs. (1) and (2)] 



s 


CO 


Cl 


C2 


C3 


C 4 


C6 





1 


20.62 


27.37 


164.19 


106.0 


1175.6 


1 


5 


19.57 


19.72 


212.0 


495.5 


1903.7 


2 


12 


12.53 


13.02 


207.8 


907.9 


1532.6 


3 


21 


3.393 


8.177 


203.6 


1313.8 


1140.1 



at each Thomson depth t arc applied itcratively until 
AT(t)/T(t) < 1%, deviations from radiative equilibrium 
are < 1%, and deviations from constant flux are < 2% (see 
HL01 for details of our numerical method). 

For the high magnetic field models (Section 4.2), we in- 
clude the effect of vacuum polarization. As shown previ- 
ously (Lai & Ho 2002, hereafter LH02), vacuum polariza- 
tion induces conversion between the low-opacity mode and 
the high-opacity mode, and the efficiency of mode conver- 
sion depends on the photon energy. In the absence of a 
rigorous treatment of this effect, we follow HL03 by con- 
sidering the two limiting cases: complete mode conversion 
(Section 4.2.2) and no mode conversion (Section 4.2.3). 
Our treatment of the no mode conversion case improves 
upon that given in HL03. 

A major complication when incorporating bound species 
in the atmosphere models arises from the strong coupling 
between the center-of-mass motion of the atom and the in- 
ternal atomic structure. PC have constructed thermody- 
namically consistent EOS and opacity models of partially 
ionized hydrogen plasmas in strong magnetic fields. From 
the work of PC, we obtain tables for the EOS and absorp- 
tion opacities of hydrogen. The EOS tables for a given 
magnetic field are arranged by temperature and density. 
The pressure at Thomson depth r (= — J pkq s dz, where 
Kq s = <7t/to p ) is calculated from hydrostatic equilibrium. 
With a given temperature profile, we search the EOS ta- 
bles for the nearest temperatures and pressures and do a 
weighted average to obtain the density profile. We also 
extract the atomic fraction from the EOS tables. 

The opacity tables contain the absorption opacity «£ bs 
(where a — 0, ±) as a function of temperature, density, 
and photon energy. The absorption opacity for the modes 
(j = 1 for X-mode, j = 2 for O-modc) is obtained from 



abs _ \ " 



2^ 



j 1 2 abs 



(3) 



where e J 



e{ is the z-component (along B) of the mode 
eigenvector and e± — (ej ± ie J y )/\/2 are the circular com- 
ponents. We calculate the polarization eigenvectors e J as- 
suming that the medium is fully ionized since the abun- 
dance of bound species is very low for the cases considered 
here (see Section 4); this is not strictly correct, but the 
effect of a small amount of bound species on the polar- 
ization of the medium is not known at present (see Bu- 
lik & Pavlov 1996 for the case where the atmosphere is 
completely neutral). 1 For a given layer in the atmosphere 
model, which specifies a particular temperature and den- 
sity, the opacity k^> s (E) at a given energy E is obtained 



from the geometric average of the nearest two or three 
energies of the table and used to compute Kj hs (E). The 
scattering opacities are computed by the same method as 
in HL01 and HL03. 



4. RESULTS 

4.1. Low-Field (B < B Q ) Models 

We consider three H atmosphere models in this sub- 
section: B 12 = 1, T off , 6 = Xcff/(10 6 K) = 0.5; B 12 = 
l,T e ff t Q = 1; B\2 = 10, Toff. 6 = 1. Figure 1 shows the tem- 
perature and atomic fraction profiles as a function of den- 
sity. The atomic fraction is the number of H atoms with 
non-destroyed energy levels divided by the total number of 
protons (see PC). The dependence of the atomic fraction 
on temperature and magnetic field is evident; in particu- 
lar, lower temperatures or higher magnetic fields increase 
the abundance of bound species, with a maximum of a few 
percent for the models shown here. 

Figure 2 shows the angle-averaged X-mode absorption 
opacity K^ s , defined by 



K X 



1 

47T 



dk Kf S (k). 



(4) 



as a function of energy for various temperatures and den- 
sities typical of the atmosphere models considered here. 
Even though the bound species abundance is very low, 
Fig. 2 shows that they can still have a very strong influ- 
ence on the opacities. Besides the proton cyclotron line at 
Eb p = 6.3i?i2 eV, there are several strong features due to 
bound-bound and bound-free transitions. These are the 
s = to s = 1 transition at E = 51 eV for B12 = 1 and 
at 0.14 keV for B12 = 10, the s = to s = 2 transition at 
E = 75 eV for B 12 = 1 and at 0.23 keV for B 12 = 10, and 
the bound-free transition at E = 161 eV for Bi 2 = 1 and 
at 0.31 keV for B 12 = 10. Because of magnetic broadening 
(see Section 2), these features resemble bumps rather than 
ordinary spectral lines. 

Figure 3 shows the spectra of the three atmosphere mod- 
els. For comparison, we also show the spectra of the fully 
ionized models at the same B and T e ff, along with the 
non-magnetic fully ionized models at the same T e g and the 
blackbody spectrum with T = T e g. The proton cyclotron 
line at Eb p — 6.3 Bi 2 eV is clear in both the fully ion- 
ized and partially ionized spectra. 2 At B < Bq, the effect 
of vacuum polarization on the atmosphere structure and 
spectra is weak; therefore there is no suppression of the 
proton cyclotron line (see Section 4.2 and HL03), and the 
proton cyclotron line can be very broad, especially when 



1 For all the models considered in this paper, the neutral fraction is less than a few percent. 

2 The width of the proton cyclotron line in all the partially ionized models is due to the finite energy grid of the models and not an indication 



of the true line width; the true width is narrower. 
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Fig. 1. — Temperature (upper panel) and atomic fraction (lower panel) as functions of density for various partially ionized H atmosphere 
models discussed in Section 4.1: magnetic field B = 10 12 G and effective temperature T e fp = 5 X 10 5 K (solid lines), B = 10 12 G and 
T e ff = 10 6 K (dotted lines), and B = 10 13 G and T e ff = 10 6 K (dashed lines). Atomic fraction is the number of H atoms with non-destroyed 
energy levels divided by the total number of protons. 




E (keV) 



Fig. 2. — Angle-averaged X-mode absorption opacity as a function of energy for temperatures and densities representative of the 

atmosphere models discussed in Section 4.1. The upper panel shows K^ B appropriate for the model with magnetic field B = 10 12 G and 
effective temperature T eS = 5 X 10 5 K — T = 4.2 X 10 5 K and p = 0.093 g cm" 3 (dotted line), T = 2.4 X 10 s K and p = 4.0 X 10~ 3 g cm" 3 
(dashed line), T = 2.0 X 10 5 K and p = 1.9 X 10 — 4 g cm~ 3 (solid line). The middle panel shows i\j^ bs appropriate for the model with 
B = 10 12 G and T cB = 10 6 K — T = 7.4 X 10 5 K and p = 0.16 g cm" 3 (dotted line), T = 6.3 X 10 5 K and p = 0.061 g cm~ 3 (dashed 
line), T = 4.3 X 10 5 K and p = 8.7 X 10~ 3 g cm -3 (solid line). The lower panel shows K J bs appropriate for the model with B = 10 13 G and 
T cff = 10 6 K — T = 8.4 X 10 5 K and p = 1.0 g cm~ 3 (dotted line), T = 6.2 X 10 5 K and p = 0.061 g cm~ 3 (dashed line), T = 3.5 X 10 5 K 
and p = 1.1 X 10 — 3 g cm -3 (solid line). The opacity curves in each panel are in descending order from hotter and higher density to cooler 
and lower density, and the upper K^ a (dotted) curves have been multiplied by 10 while the lower K^? s (solid) curves have been divided by 
10 for clarity. See text for description of the opacity features. 
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Fig. 3. — Spectra of hydrogen atmospheres with B = 10 12 G and T e g = 5 X 10 5 K (upper panel), B = 10 12 G and T e g = 10 6 K (middle 
panel), and B = 10 13 G and T cf f = 10 6 K (lower panel). The solid line is for a partially ionized atmosphere, the dot-dashed line is for a fully 
ionized atmosphere, the dashed line is for a fully ionized non- magnetic atmosphere, and the dotted line is for a blackbody with T = T e g . The 
proton cyclotron line in the partially ionized spectra is not resolved (see footnote 2). 



E Bp > 3kT cS (HL01). A comparison of the B = 10 12 G 
models shows that spectral features due to neutral species 
in the atmosphere are weaker in the higher temperature 
model. This is because the higher temperatures decreases 
the abundance of neutral species (see Fig. 1) and thus re- 
duces the effect of the neutral species on the continuum 
opacity. 

It is evident from Fig. 3 that spectral features due to 
bound species at these field strengths lie within the ex- 
treme UV to very soft X-ray (50 eV < E < 200 eV) regime 
and therefore are difficult to observe due to interstellar 
absorption. However, the effect of bound species on the 
temperature profile of the atmosphere and the continuum 
flux is significant. In particular, the optical flux is higher 
for the partially ionized atmospheres compared to the fully 
ionized atmospheres. On the other hand, the partially and 
fully ionized models both yield very similar (neglecting the 
spectral features) hard X-ray flux for a given effective tem- 
perature, with the partially ionized model fluxes slightly 
lower than those from the fully ionized models; thus fit- 
ting the observed hard X-ray flux with fully ionized models 
would yield fairly accurate NS temperatures. 

4.2. Magnetar (B > Bq) Models 

When studying magnetar atmosphere models, it is im- 
portant to include the effect of vacuum polarization. How- 
ever, to clearly see the role of bound species, we first con- 
sider models in which vacuum polarization is neglected. 

4.2.1. No Vacuum Polarization 

Figure 4 shows the temperature and atomic frac- 
tion profiles for the H atmosphere model with _B 14 = 
£>/(10 14 G) = 5 and T e ff,6 = 5. Here we use an exten- 
sion of the PC EOS and opacity model to B > 10 13 ' 5 G. 



In this case, instead of the fits of Potekhin (1998) employed 
in PC (some of which would fail at B > 10 13 ' 5 G), we per- 
form numerical calculations of the atomic energies, sizes, 
and oscillator strengths following Potekhin (1994). (Some 
results of this extension of the PC model have been re- 
ported in Chabrier, Douchin, & Potekhin 2002.) Though 
the magnetic field increases the binding energy of the hy- 
drogen atom to > 0.5 keV, the atomic abundance remains 
< 1% in the atmospheric layers. 

The upper panel of Figure 5 shows the angle-averaged 
X-mode opacity K^? s as a function of energy for various 
temperatures and densities typical of the high magnetic 
field model considered here. Besides the proton cyclotron 
line at Eb p = 0.63 -B14 keV, the sharp change in opacity 
at 0.76 keV for B14 = 5 is due to photoionization. The 
features at E > Eb p are the partial photoionization edges 
for transitions to the continuum states with s > 0, which 
are merged with autoionization resonances just below the 
corresponding partial thresholds and extended redwards 
due to magnetic broadening (see Potekhin & Pavlov 1997). 
The bound-bound transitions are not seen in Fig. 5 be- 
cause the bound states with s > merge with the contin- 
uum at such field strengths because of the sHuj Bp-term in 
E s . 

Figure 6 shows the spectra of the models with B = 
5 x 10 14 G and T cS = 5 x 10 6 K. The upper panel is 
a comparison of partially and fully ionized models with 
no vacuum polarization. Broad absorption features at 
E ~ 0.76 keV and 4 keV are due to bound-free transi- 
tions to different continuum states (the latter feature is 
from the blending of two magnetically broadened partial 
photoionization edges at 3.4 keV and 6.5 keV; see Fig. 5), 
while the continuum flux between the models is not sig- 
nificantly different. Also, as noted in HL01 and HL03 (see 
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Fig. 4. — Temperature (upper panel) and atomic fraction (lower panel) as functions of density for the partially ionized H atmosphere 
models discussed in Section 4.2: B = 5 X 10 14 G and T off = 5 X 10 6 K with vacuum polarization (solid lines) and without vacuum polarization 
(dashed lines). Atomic fraction is the number of H atoms with non-destroyed energy levels divided by the total number of protons. 




E (keV) 

Fig. 5. — Angle- averaged absorption opacities as functions of energy for temperatures and densities representative of the B = 5 X 10 14 G 
and T e ff = 5 X 10 6 K atmosphere model including vacuum polarization discussed in Section 4.2. The upper panel shows the X-mode opacity 
K^ B at T = 2.9 X 10 6 K and p = 0.13 g cm~ 3 (dotted line), T = 2.4 X 10 6 K and p = 0.031 g cm" 3 (dashed line), T = 1.8 X 10 6 K and 
p = 4.0 X 10 — 3 g cm~ 3 (solid line). The opacity curves in the upper panel are in descending order from hotter and more dense to cooler and 
lower density, and the upper K^ B (dotted) curves have been multiplied by 10 while the lower K^ 11 (solid) curves have been divided by 10 
for clarity. The lower panel shows the plus-mode and minus-mode opacities K^ ba (J = ±) and T and p corresponding to the dotted curve in 
the upper panel; using K^ s and K^° s is appropriate for models in which complete mode conversion is assumed (see Section 4.2.2). See text 
for description of the opacity features. 
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Fig. 6. — Spectra of hydrogen atmospheres with B = 5 X 10 G and T e ff = 5 X 10 K. The upper panel shows the partially ionized 
(solid line) and fully ionized (dashed line) atmospheres with no vacuum polarization. The lower panel shows the partially ionized (solid line) 
and fully ionized (dashed line) atmospheres with vacuum polarization and assuming complete mode conversion. The dotted lines are for a 
blackbody with T = 5 X 10 6 K. The proton cyclotron line in the partially ionized spectra is not resolved (see footnote 2). 



also Zane et al. 2001), the proton cyclotron line is very 
broad when vacuum polarization is neglected. 

4.2.2. Including Vacuum Polarization: Complete Mode 

Conversion 

The main effect of vacuum polarization on radiative 
transfer is the "vacuum resonance," which occurs when 
the effects of the vacuum and plasma on the polarization 
modes cancel each other. For a photon of energy E{= hui), 
the vacuum resonance occurs at density 

p v « 0.96 Y^ElBljiB)- 2 g cm" 3 , (5) 

where E\ = E/ (1 keV) and f(B) is a slowly- varying func- 
tion of B of order unity. The two photon modes are almost 
linearly polarized away from the resonance, while they 
both become circularly polarized near py. LH02 showed 
that a photon propagating in the atmospheric plasma can 
adiabatically convert from the low-opacity X-mode to the 
high-opacity O-mode (and vice versa) as it crosses py. For 
this conversion to be effective, the adiabatic condition 

#ad > 1-49 [/(S)tan0 B |l - Ui |] 2/3 (5 cm/H p ) 1/3 keV, 

where 9b is the angle between k and B, u\ = iv Bp /uj 2 , and 
Hp = \dz/d\np\ is the density scale height along the ray, 
must be satisfied. Clearly, resonant mode conversion can 
significantly affect radiative transfer (LH02; HL03; Lai & 
Ho 2003). ' 

To rigorously treat the radiative transfer across the vac- 
uum resonance, one must go beyond the modal descrip- 
tion of the radiation field by solving the transfer equation 
in terms of the photon intensity matrix (Lai & Ho 2003). 
In the present paper, as in HL03, we consider two lim- 
iting cases. In the "complete mode conversion" case, we 



assume photons of all energies change their mode charac- 
teristics (from X to O-mode or vice versa) when crossing 
the vacuum resonance. This can be achieved by using the 
plus- mode and minus-mode basis (j = ±) in the radia- 
tive transfer equation (see HL03 for details). At a given 
density p, the characteristics of the ±-modes switch at the 
energy 

Ey « 1.02 Y^ 2 p\ /2 B^f(B) keV, (7) 

where p\ = p/(l g cm -3 ). 

Figure 4 shows the temperature and atomic fraction pro- 
files for the partially ionized model which includes vacuum 
polarization in the complete mode conversion limit. As 
noted in HL03, vacuum polarization increases the tem- 
perature at each depth relative to the case without vac- 
uum polarization. The higher temperatures result in lower 
atomic fractions for a given magnetic field. 

The lower panel of Fig. 5 shows the angle-averaged 
opacity Kj hs for j = ± as a function of energy at p = 
0.13 g cm -3 and T — 2.9x 10 6 K, which are characteristic 
of the B = 5 x 10 14 G atmosphere. The distinction be- 
tween the plus and minus-modes (j = ±) and the X and 
O-modes is discussed in HL03 (in particular, see Fig. 4 of 
that work). At this particular temperature and density, 
the plus-mode (minus-mode) is in the high (low) opacity 
state for E < Ey but manifests as the low (high) opacity 
state for E > Ey. 

Fig. 6 compares the spectra of different models with 
and without vacuum polarization. First, we see that, with 
vacuum polarization, the continuum flux at high ener- 
gies for the partially ionized model is very similar to the 
fully ionized model [for the parameters shown in Fig. 6, 
F v (partial)/^ (full) < 1.5 at E > 0.1 keV]. Second, vac- 
uum polarization softens the high energy tail of the spec- 
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tram in both the partially ionized and fully ionized cases, 
although the emission is still harder than the blackbody 
spectrum because of the non-grey opacities. The softer 
spectrum is due to the density dependence of the vac- 
uum resonance energy Ey and the large density gradient 
present in the atmosphere; the change in opacity at the 
vacuum resonance causes the high energy photons emerg- 
ing from the atmosphere to be produced from a shallower, 
cooler layer in the atmosphere than in the case when vac- 
uum polarization is neglected (see Fig. 2 of Ho et al. 2003). 
Third, as discussed in HL03, this vacuum-induced depres- 
sion of continuum flux strongly suppresses the proton cy- 
clotron line by making the photospheres inside and outside 
the line more similar (see also footnote 2). We now see that 
vacuum polarization also suppresses spectral features due 
to the bound species of hydrogen. The reduced width of 
the proton cyclotron line and spectral features associated 
with bound transitions makes these features difficult to 
observe with current X-ray detectors (see Section 6). 

4.2.3. Including Vacuum Polarization: No Mode 
Conversion 

The opposite limit for treating the vacuum polarization 
effect is to assume no mode conversion across py. This is 
achieved by using the X-mode and O-mode as the mode 
basis in the transfer equation (see HL03). For this case, 
the X-mode opacity has a much lower opacity than the O- 
mode opacity except at the resonance peak E = Ey (where 
Kx = «o). As discussed in LH02, even without mode con- 
version, the optical depth across this opacity feature can 
be significant; therefore, the X-mode radiative transfer is 
affected. 

We present here improvements made to our models with 
no mode conversion. The atmosphere models contain en- 
ergy grid points E n and density grid points pi. As de- 
scribed in HL03, we use an adaptive energy grid to ac- 
count for the narrow, density-dependent opacity feature 
at Ey, i.e., an equal number of density and energy grid 
points is used with every energy grid point being placed 
at E n = Ey(pi). We noted in HL03, however, that this 
implementation overestimates the strength of the opacity 
feature on the atmosphere structure and spectra. This is 
illustrated in Figure 7: The area under the short-dashed 
line overestimates the area under the solid line, which rep- 
resents the true opacity. 

To rectify this problem, we introduce an effective opac- 
ity at the resonance, R(Ey), given by 



R{E n = Ey) 



n+l 



E, 



n-l 



k(E) dE 



1 



K(E n +i) (E n +i — E n ) 



(8) 



where the integral is calculated using an auxiliary energy 
grid that has much finer resolution. Using R(Ey) preserves 
the equivalent width of the resonance feature: In Fig. 7, 
the area under the long-dashed line is equal to the area 
under the solid line. As would be expected, atmosphere 
models based on this improved method but with different 
numbers of energy grid points yield very similar results. 



oc 



The spectra for fully ionized hydrogen atmosphere mod- 
els with B = 5 x 10 14 G and T cB = 5 x 10 6 K is shown in 
Figure 8. Note that the width of the proton cyclotron line 
in the (no mode conversion) adaptive energy grid models 
is not the true width of the line. The broad depression in 
the spectrum due to the density dependence of the vac- 
uum resonance feature is shifted to higher energies since 
the old calculation (see HL03) overestimated the effect of 
the resonance feature. The two limiting cases of complete 
mode conversion and no mode conversion should approxi- 
mately bracket the true spectrum. We also note that the 
no mode conversion flux spectrum should be suppressed 
at low energies due to the dense plasma effect described 
in Section 5, and thus the two limiting cases should have 
more similar fluxes at lower energies than is shown in Fig- 
ure 8. 

4.2.4. Proper Treatment of Free- Free Opacity 

As discussed in PC, the traditional formula (see Pavlov 
et al. 1995) for the opacity due to free- free absorption of 
both electrons and ions given by 

ff 1 f m e \ 2 1 /. x 

n a oc 2 + 1 I 2 (incorrect), 

(uj + aujBe) \ m pJ {ijJ - awBp) 

(9) 

where a = 0,±1 and Hlubc = 11.6 B12 keV is the elec- 
tron cyclotron energy, misses an interference effect due 
to electron- ion collisions. The correct free- free absorption 
opacity is given by (see, e.g., PC for more discussion) 
to 2 

-j— -2 (correct). (10) 

[U! + aUJ B e ) (W - CtUBp ) 

When lu lob p , equations (9) and (10) are very nearly the 
same. However, when lu < wb p , the correct opacity is sup- 
pressed by a factor ~ uo 2 /uo 2 Bp compared to the incorrect 
one. Thus we expect that, in the atmosphere models us- 
ing the correct opacity given by equation (10), the photons 
with E < Eb p will be generated from deeper, hotter layers 
as opposed to models which use the incorrect opacity given 
by equation (9). All previous atmosphere models adopted 
equation (9) and are therefore problematic if photons with 
E < Eb p are important (see, however, Lloyd 2003). 

Figure 9 compares atmosphere spectra generated using 
the correct expression for the free-free absorption opac- 
ity [eq. (10)] and using the incorrect expression [eq. (9)]. 
The model shown here is a more extreme example since 
Eb-p <] kT e ff (in order to maintain constant total flux, 
a large redistribution of the specific energy flux occurs). 
When Eb p <C kT c s, such as for the low magnetic field 
cases considered in Section 4.1, the difference between us- 
ing the correct expression versus the incorrect one is small 
(except possibly in the region to < u>b p ] but this region of 
the spectrum carries negligible flux). All models consid- 
ered in this paper use the correct expression. 

5. DENSE PLASMA EFFECT 

Here we comment on a main uncertainty associated with 
atmosphere models of strongly magnetized NSs. Although 
we have included the medium effect on the polarization 
modes, at sufficiently high densities or low photon ener- 
gies, the refractive index of the medium deviates appre- 
ciably from unity. In particular, at the photon decou- 
pling layer in the atmosphere (where the optical depth 
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Fig. 7. — Schematic picture of the radiative opacity as a function of energy at the density grid point pi. The vacuum resonance feature 
occurs at higher energies at larger densities (depths) [see eq. (7)]. In the adaptive energy grid method of HL03 and Section 4.2.3, the energy 
grid points E„_i, E n , E n +i are chosen to be equal to the vacuum resonance energy at one of the density grid points, i.e., E n = Ey(pi). The 
short-dashed line indicates the apparent opacity obtained by interpolating energy grid points E n -i, E n , -En+i, while the long-dashed line 
indicates the apparent opacity obtained by the method described in Section 4.2.3. 
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Fig. 8. — Spectra of fully ionized hydrogen atmospheres with B = 5 X 10 14 G and T c ff = 5 X 10 6 K. The solid line and dashed line are 
the atmospheres with vacuum polarization but no mode conversion and complete mode conversion, respectively, the dot-dashed line is the 
atmosphere without vacuum polarization, and the dotted line is for a blackbody with T = 5 X 10 6 K. 
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r„ i=a 1), the plasma frequency ui pe may exceed the pho- 
ton frequency u>. For example, in the B = 10 12 G and 
T e s = 5 x 10 5 K model shown in Section 4.1, the X-mode 
photons (which are the main carriers of the energy flux) 
with E < 10 eV decouple from the atmosphere layer where 
E < huj pe = 28.71 (p/1 g cm" 3 ) 1 / 2 cV (see lower panel of 
Fig. 10). Clearly we expect that the actual emission at 
these energies is lower than what is shown in our model 
spectrum (see upper panel of Fig. 3). 

A quantitative study of this dense plasma effect is be- 
yond the scope of this paper. However, we can obtain an 
estimate of its effect by a simple procedure. In the atmo- 
sphere model, we evaluate the parameter v c (— lo^/lo 2 ); 
whenever v c > 1, we artificially increase all the opacities 
at that energy and depth by a factor 7 P (S> 1). This sim- 
ulates the suppression of radiation below the plasma fre- 
quency cutoff. The resulting spectrum is shown in Fig. 10, 
where we have chosen j p = 10 2 and 10 4 and labeled the 
results as "dense." The result of Section 4.1, where we 
"blindly" carry out the calculations even when v e > 1, is 
labeled as "normal" in Fig. 10. We see that the spectrum is 
largely unaffected since v e < 1 at most energies. However, 
we see the flux at low energies (E < 5 eV) is suppressed 

[with i^^/^dense) ^ 1A ^ ^ ^ for ^ = ^ ^ 
F („orm) /i? (dcnsc) ^ ^ ^ ^ ^ ^ ^ = ^ ^ ^ 

increase in flux at ~7-20 eV is due to flux redistribution 
to maintain constant total flux or radiative equilibrium. 

For the magnetar models, the strong suppression of the 
X-mode opacities by the magnetic field causes the pho- 
tons to be generated from deeper, denser atmosphere lay- 
ers. Therefore, the dense plasma effect may be even more 
important. (However, this can be offset by vacuum polar- 
ization which makes the photosphere shallower; see Sec- 
tion 4.2.) Clearly this dense plasma effect needs to be un- 



derstood better before accurate comparisons can be made 
with, for example, optical and infrared observations. 

6. DISCUSSION 

We have presented results of the first magnetic hydro- 
gen atmosphere models of isolated neutron stars based on 
the best current equation of state and opacities of magne- 
tized, partially ionized H plasmas (PC). For models with 
"ordinary" magnetic field strengths (10 12 G< B < Bq = 
4.4 x 10 13 G), the spectral lines associated with bound 
species lie in the extreme UV to very soft X-ray energy 
bands and are difficult to observe due to interstellar ab- 
sorption at these energies (though a broad absorption fea- 
ture at ~ 0.2 keV was seen in RBS 1223; Haberl et al. 2003; 
see also Ho & Lai 2003b). However, the opacities are suffi- 
ciently different from the fully ionized opacities that they 
can change the atmosphere structure and continuum flux. 

We find that, for the higher magnetic field models 
(B > -Bq), vacuum polarization suppresses both the pro- 
ton cyclotron line and spectral lines due to the bound 
species. As a result, the thermal spectra are almost fea- 
tureless and blackbody-like. Recent observations of several 
AXPs with Chandra and XMM-Newton X-ray telescopes 
failed to resolve any significant line features in the spectra 
(e.g., Patel et al. 2001, 2003; Juett et al. 2002; Tiengo et 
al. 2002; Morii et al. 2003); thus this could possibly be 
indicating the presence of vacuum polarization effects in 
these magnetar candidates. 

Our results may also be useful in understanding spec- 
tral observations of dim isolated neutron stars, such as 
RX J1856. 5-3754 and RX J0720.4-3125 (see, e.g., Pavlov 
et al. 2002; Kaplan et al. 2003). Since kT cS ~ tens of eV 
for these neutron stars, the presence of neutral atoms must 
be accounted for in the atmosphere models if B > 10 12 G. 
We find that one-temperature partially ionized hydro- 
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Fig. 10. — Spectra (upper panel) and density at the photosphere (lower panel) for partially ionized hydrogen atmospheres with B = 10 12 G 
and T cf f = 5 X 10 5 K. In the upper panel, the solid line shows the spectrum ("normal") when the dense plasma effect discussed in Section 5 
is neglected, while the dashed line and dotted line show the spectra ("dense") when the dense plasma effect is approximately accounted for 
using 7 P = 10 2 and 10 4 , respectively. The lower panel shows the density at the X-mode photon decoupling layer (at which t„ tts 1; solid line) 
and the O-modc photon decoupling layer (dashed line). The dotted line denotes the density at which the photon frequency equals the plasma 

1 /2 

frequency, i.e., E = Ijuipe = 28.71 p-j/ eV. 



gen atmosphere models improve the fitting of the X-ray 
data (compared to fully ionized models), although they 
still do not provide a satisfactory fit to the combined X- 
ray and optical spectra. One possible solution is "thin" 
atmosphere models in which a gaseous hydrogen atmo- 
sphere lies on top of a condensed matter (see Motch et 
al. 2003). We find that the spectra of the neutron stars 
RX J1856.5 - 3754 and RX J0720.4 - 3125 can be crudely 
fit by our B = 10 12 G and T off ~ 5 x 10 5 K partially ionized 
hydrogen "thin" atmosphere model (Ho et al., in prepa- 
ration). The resulting optical spectrum of the model only 
overpredicts the optical data by < 1.5. Furthermore, when 
the dense plasma effect discussed in Section 5 is taken into 
account, the suppression of low energy flux may lead to an 
even better fit of the optical spectrum. 

Chandra and XMM-Newton observations revealed ab- 
sorption features at 0.7 and 1.4 keV in the spectra of the 
neutron star IE 1207.4 - 5209 (Mereghetti et al. 2002b; 
Sanwal et al. 2002). Sanwal et al. (2002) contend the fea- 
tures are due to atomic transitions of singly-ionized helium 
in the neutron star atmosphere with a superstrong mag- 
netic field, B w 1.5 x 10 14 G. Hailey & Mori (2002) and 
Mori & Hailey (2003) attribute the features to transitions 
of mid-Z oxygen or neon in a ~ 10 12 G magnetic field, 
while Mereghetti et al. (2002b) interpret the features to be 
the result of high-Z metals at - 10 12 G. If IE 1207.4-5209 
is indeed a neutron star with superstrong magnetic field 
(B > Bq) and the properties of the helium atom are simi- 
lar to those of the hydrogen atom at these field strengths, 
then our results show that any spectral features would not 
be resolved due to their suppression by vacuum polariza- 
tion. However, we cannot be more definitive because of the 



unknown properties of helium at these field strengths and, 
as discussed in Section 3, we assumed that the presence 
of a very small amount of neutral atoms has no effect on 
the polarization properties of the medium. Furthermore, 
recent observations have revealed an additional feature at 
2.1 keV and possibly another at 2.8 keV thus hinting at an 
electron cyclotron origin (with B = 8 x 10 10 G) for these 
features (Bignami et al. 2003). 

Several caveats/uncertainties are worth mentioning. (1) 
The radiative transfer formalism adopted in our work relics 
on the transfer of two photon polarization modes. This is 
inadequate for treating the vacuum polarization-induced 
resonant mode conversion effect, especially because the ef- 
fectiveness of mode conversion depends on photon energy 
(see LH02; HL03). To properly account for the mode con- 
version effect, as well as mode collapse and the breakdown 
of Faraday depolarization, one must go beyond the modal 
description of the radiation field by formulating and solv- 
ing the transfer equation in terms of the photon intensity 
matrix (or Stokes parameters) and including the effect of 
a nontrivial refractive index (see Lai & Ho 2003). (2) As 
discussed in Section 5, the dense plasma effect on radiative 
transfer can be important, especially for low energy pho- 
tons and high magnetic fields. These issues require further 
study. 
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